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Abstract 

^ The magnetic field in many astrophysical plasmas, for example in the solar 

corona, is known to have a highly complex - and clearly three-dimensional - 
structure. Turbulent plasma motions in high-/? regions where field lines are 
anchored, such as the solar interior, can store large amounts of energy in the 
P^h magnetic field. This energy can only be released when magnetic reconnection 

i— h occurs. Reconnection may only occur in locations where huge gradients of 

the magnetic field develop, and one candidate for such locations are magnetic 
p/j null points, known to be abundant for example in the solar atmosphere. Re- 

C/} connection leads to changes in the topology of the magnetic field, and energy 

43 released as heat, kinetic energy and acceleration of particles. Thus reconnec- 

tion is responsible for many dynamic processes, for instance flares and jets. 
The aim of this paper is to investigate the properties of magnetic reconnec- 
tion at a 3D null point, with respect to their dependence on the symmetry of 
the magnetic field around the null. In particular we examine the rate of flux 
transport across the null point with symmetric/asymmetric diffusion regions, 
as well as how the current sheet forms in time, and its properties. Mathemat- 
ical modelling and finite difference resistive MHD simulations. It is found 
that the basic structure of the mode of magnetic reconnection considered is 
unaffected by varying the magnetic field symmetry, that is, the plasma flow 
is found cross both the spine and fan of the null. However, the peak intensity 
and dimensions of the current sheet are dependent on the symmetry/ asym- 
metry of the field lines. As a result, the reconnection rate is also found to 
be strongly dependent on the field asymmetry. The symmetry /asymmetry 
of the magnetic field in the vicinity of a magnetic null can have a profound 
effect on the geometry of any associated reconnection region, and the rate at 



(N 
> 



X 



which the reconnection process proceeds. 



1 Intoduction 



Magnetic reconnection is the breaking and topological or geometrical rear- 
rangement of the magnetic field lines in a plasma. The magnetic field plays 
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a fundamental role in many of the phenomena that occur in the plasma. It is 
not surprising that three-dimensional (3D) magnetic fields are more complex 
than two-dimensional ones. It is known from observations that magnetic 
reconnection occurs in abundance in astrophysical plasmas. However, due 
to the very low plasma resistivity, reconnection may only occur where very 
intense currents ('current sheets') develop. One of the most fundamental 
questions that must be answered to determine the locations and mechanisms 
of energy release in astrophysical plasmas is therefore: where may such cur- 
rents develop? 

It is now becoming clear that the magnetic field in the solar atmosphere 
has a highly complex structure. Two major candidates that have been pro- 
posed as sites of current sheet formation in such a complex magnetic field 



are 3D nulls points, and associated separator field lines (Lau and Finn 



Klapper et al. 1996 Priest and Titov 1996 Longcope and Cowley 



Galsgaard and Nordlund 1997 Pontin and Craig 2005 Longcope 





1990 




1996 




1996 



2001) - field lines that link two nulls. We focus here on reconnection at iso- 



lated 3D nulls. Indications are that an abundance of 3D nulls is present in the 



solar corona (Regnier et al. , 2008 Longcope and Parnell, 2009), which have 



been suggested as likely sites for coronal heating (e.g. Priest et al. 2005). 



Moreover, recent observations suggest that reconnection at such nulls may 



play an important role in jets (|Pariat et al. 
flares 



(e.g. 



Luoni et al. 



e.g. Ugarte-Urra et al. 



2007 



2007 



Masson et al. 



2009; Torok et al., 2009) solar 



2009) and coronal mass ejections 



Barnes, 2007). Furthermore, recently the first 



in-situ observations have been made by the Cluster spacecraft of single and 



multiple 3D magnetic nulls in the Earth's magnetotail (e.g. Xiao et al. , 2006). 
The observations further suggest that these nulls may be playing an impor- 
tant role in the reconnection process occurring in the magnetotail. Though 
magnetic field measurements in other astrophysical objects further afield are 
difficult, it is almost certain that similar reconnection processes at 3D nulls 
also occur there. 

To find the local magnetic structure about a null point, we consider the 
magnetic field in the vicinity of null point where the field vanishes (B = 0). 
If the null point is taken to be situated at the origin and, in addition, we 
assume we are sufficiently close to the null, then the magnetic field may be 
expressed as 

B = Mr (1) 

where Ai is a matrix with the elements of the Jacobian of B and r is the 
position vector (x, y, z) T . The eigenvalues of Ai sum to zero since V ■ B = 0. 
We consider the situation where all the eigenvalues are real. Since they sum 
to zero there is always one eigenvalue of opposite sign to the other two. 
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The two eigenvectors corresponding to the eigenvalues with same-sign real 
part define the "fan surface" of the null. The third eigenvector defines the 
orientation of the "spine line". For more details, see e.g. Fukao et al. (1975); 



Lau and Finn (1990); Parnell et al. (1996). 



Unlike in two dimensions, reconnection can occur in 3D either at a null 



point or in the absence of nulls (Schindler et al. 1988 Priest and Forbes 



2000; Demoulin 2006). What's more, the nature of reconnection in 3D 



has been shown to be fundamentally different from 2D reconnection (Priest 



et al. , 2003). The nature of magnetic reconnection in the absence of a three- 



dimensional null point has been discussed by Hesse (1991) and Hornig and 



Priest (2003). The kinematics of steady reconnection at three dimensional 
null points have been studied by Priest and Titov (1996) when rj — 0. Later, 



Pontin et al. (2004, 2005) improved this model by adding a finite resistivity, 
localised around the null point. Two distinct cases were considered, in which 
the current (J) was directed parallel to first the spine and second the fan 
plane of the null. The structures of the two solutions were found to differ 
greatly, and as a result, the reconnection rate, calculated by integrating the 
E\\ along field lines, represents very different behaviors of the flux for the 
two cases. In the first case, in which J was directed parallel to the spine, a 
type of rotational flux mis-matching was found, with no flow being present 
across either the spine or the fan of the null point. On the other hand, when 
J was directed parallel to the fan surface, it was found that magnetic flux 
is transported through the spine line and the fan plane, in a process much 
more conceptually similar to the 2D case. In this case it can be shown that 
the reconnection rate gives a measure of the rate of flux transport across the 
separatrix surface of the null (Pontin et al. 2005). The case in which J 



is 



parallel to the spine corresponds to one pair of complex conjugate eigneval- 
ues, whereas when J is parallel to the fan the eignevalues are all real. In each 
of these investigations only the azimuthually symmetric case was considered, 
that is the case in which the magnetic field in the fan plane is isotropic. 
In this paper we focus on the case where J is parallel to the fan surface 
(real eigenvalues), and for the first time consider magnetic reconnection at 
a generic non-symmetric magnetic null point, i.e. a null for which the fan 
eigenvalues are not equal. The different modes of reconnection that occur 
in practice in a plasma (when the full set of MHD equations are considered) 
have recently been classified by Priest and Pontin (2009). In terms of the 
framework they have set up, the mode of reconnection considered here is 
termed spine-fan reconnection. In a future paper we will go on to generalise 
the complex conjugate eigenvalue case. 

In sections [2] and [3j we describe a kinematic model for reconnection at a 



non-symmetric null point, comparing our results with those of Pontin et al. 
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(2005). In section 111 we describe the results of a related resistive magneto- 



hydrodynamic (MHD) numerical simulation, and in section [5] we present our 
conclusions. 



2 Kinematic solution — method 



2.1 The model 

The subject of magnetic reconnection is a complex one, and its study is 
still in the early stages. Therefore, one approach that is used to try to 
understand the properties of this process is to consider a reduced set of 
the MHD equations. There are a number of analytical 3D solutions, which 



are described by Hornig and Priest (2003) and Wilmot-Smith et al. (2006 



2009), where there is no null point of the magnetic field, as well as the 



solutions in the presence of a null mentioned above (Pontin et al. , 2004, 2005 



Priest and Pontin, 2009). These solutions are kinematic reconnection, that 



is they satisfy Maxwell's equations, as well as the induction equation. This 
approach can give great insight into the topological structure of a magnetic 



reconnection process occurring at an isolated diffusion region (Schindler et al. 



1988). After investigating the properties of the solutions of this subset of the 



MHD equations, we go on in section [4] to examine which properties survive 
when the full set of resistive MHD equations is solved. 

We seek a solution to the kinematic, steady, resistive MHD equations in 
the locality of a magnetic null point. That is, we solve 



E + v x B = r] J 
V x E = 
V x B = fi 3 
VB = 



(2) 
(3) 
(4) 
(5) 



As discussed above, here we consider a null point with current directed par- 
allel to the fan plane. We choose the magnetic field to be 



B 



Bn 



L p + 1 



(x,py-jz, -{p+ l)z) 



(6) 



where p is a parameter (here we restrict ourselves to the case p > 0). This 
generalises the previous work by Pontin et al. ( |2005 ), who considered only 
the case where the field in the fan plane (z = 0) is azimuthally symmetric, 
corresponding to p = 1. For convenience we will write 2B Q /L(p + 1) = B' . 
The current lies in the x-direction, and is given by J = (B' Q / p,o)(j, 0, 0) from 
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Eq. Q. Examining the matrix M. (see Eq. [T]), the eigenvalues of the null 
point are found to be 



Ai = ^, X 2 =pB' , \ 3 = -(p+l)B' 
with corresponding eigenvectors 

k 1 = (1,0,0), k 2 = (0,1,0), k3=fo,i,^±i) 

It is clear from the above that the fan plane is defined by ki and k 2 (since 
p > 0). The fan plane of this magnetic null point is coincident with the 
plane z = while the spine is not perpendicular to this, but rather lies along 
x = 0,y = jz/(2p + 1) (see figure [TV. 

For the chosen magnetic field (ph, closed-form expressions for the equa- 
tions of magnetic field lines can be found, by solving 

^ = B(X(.)), (7) 
where the parameter s runs along field lines, to give 

x = Xo e B '° s (8) 

y=(y - -^-) e^ s + J\ e - B ^+ l > (9) 
y \ yo 2p + 1 J 2p+l v ' 

z = Zo e- B ^ +1 >. (10) 
The inverse of equations ( 8f9"fl~0 ) are 



XQ = xe -B' s (11) 

y =(y- -^—] e- pB '° s + J^— e B '^ a (12) 
yo \ y 2p + 1 J 2p + l K J 

ZQ = ze B '^ +1 > (13) 

which describes the equations of the magnetic field lines in terms of some 
initial coordinates X = (xq, yo, Zq). 

We proceed to solve (|2]{5]) as follows. From equation ^ we can write, 
in general E = — V0 where is a scalar potential. Then the component 
of equation (J2| parallel to B is — (V0)|| = 77 Jj | and we can calculate 4> by 
integrating along magnetic field lines: 



J 77 J + (14) 



5 





Figure 1: The structure of the magnetic null point with j = 1 and different values of p: 
(a) p = 0.5, (b) p = 1 and (c) p = 2. 
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where <po is a constant of integration. By substituting the equations (8j9jl0) 
into the integrand of Eq. (14), we can perform this integration to obtain 



</>(X , s). One this is done, we use equations ( JIT|l2|13 ) to eliminate s, x and 
yo to obtain 0(X), treating zo as a constant (see below). The electric field 
can subsequently be found from 



E = -Vcj) 



(15) 



and we then find the plasma velocity perpendicular to the magnetic field Vj_, 
by taking the vector product of equation ^ with B to obtain 



V_L 



(E — 77J) x B 

B~ 2 



(16) 



Now, in order to investigate the properties of magnetic reconnection in 
a fully 3D system, we impose a resistivity model which ensures that the 
diffusion region is spatially localised in 3D. This is also the case relevant to 
astrophysical plasmas, which are known to be effectively ideal except in very 
small regions where energy release occurs. The diffusion region is chosen 
to be localised around the null point, in line with the results of past work 
which has shown that shearing motions tend to focus current in the vicinity 



of the null (Rickard and Titov, 1996 Pontin and Galsgaard, 2007 Pontin 



et al. 2007). Since the current is uniform in our simple model, we choose 



the resistivity to be localised, and take it to be of the form 



V = Vo\ 



(R* _ 1 

\a 2 /P l 



_2p_ 



b 2 



(z 2 )F+T 
b 2 



R 2 < a 2 , (z 2 )p+ 1 < b 



R 2 <a 2 , (z 2 )pTT <b 2 <M 







otherwise 



where R x = ^{x 2 Y + (y - jz/(2p + l)) 2 and R 2 = ^x 2 + {{y - jz/(2p + I)) 2 ) 1 /?, 
where a and b are positive constants. This is done in order to localise the 
product rjJ, and hence the diffusion region, since we have not yet discovered 
a way to proceed with our analytical method with localised J. The exact 
mathematical form for r) is not expected to affect the qualitative structure of 
the solution, and is chosen in order to render the equations tractable. The 
crucial property for the structure of the solution is the localisation of the 
diffusive term 77 J. The dependence of 77 on p is chosen differently for p > 1 
and p < 1 to ensure that i](x,y,z) is always differentiable, and is chosen 
in such a way as to maintain consistency in the dependence of the diffusion 
region size in the x-direction on p, which is shown later to be an important 
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property. r] is the value of rj at the null point, and the diffusion region is a 
tilted cylinder centered on the spine axis, extending to z — ±fe(p +1 )/ 2 P when 
p > 1 and z = ±6(p +1 )/ 2 when p < 1. The cross-section of the diffusion 
region in the z = plane is circular with radius a when p = 1, but when 



p 1, it extends to x = ±a l l p and y = ±a. In order to integrate Eq. (14), 
we must choose a surface on which to start our integration (i.e. on which 
to set s = 0) that intersects each field line once and only once, in order 
that is single-valued. We choose surfaces above and below the fan surface, 
z = ±Zq) constant. To simplify the mathematical expressions, and without 
loss of generality, we assume zq = b. Performing the calculation of 0(X) as 
described above yields two expressions for 0, for z > and z < 0. In order 
to match these two expressions at the fan plane, that is for to be smooth 
and continuous, and thus physically acceptable, we must set the value of 



at z — ±z (i.e. O in Eq. 14) to be 

2(p-l) 

... ffo H-i -1, 0<p< 1, 

00 = 01^0 \ -, 4(p-lj „ 2(p-l) 

6 P+ 1 - 1, p > 1, 




where 



(8p-l)" (4p-l) 1 

L(p+ 1) 



(19) 



0(X), E and vj_ can be obtained from ( 14 ) ,([15]) and (16), as described ear- 
lier. The mathematical expressions are too lengthy to show here but can be 
calculated using a symbolic computation package. Here we have used Maple 
v.12. 



3 Kinematic solution — analysis 
3.1 Nature of reconnection 

In order to determine the structure of the magnetic reconnection process, we 
will examine the plasma velocity perpendicular to the magnetic field (vj_). 
This velocity transports the magnetic flux outside the diffusion region. The 
flow does not cross the spine in the x-direction (v± x (0,y, z) = 0), so that 
v± x is negligible for the reconnection process. However, in the yz-plenae, the 
plasma flow crosses both the spine and the fan. Note that this is qualitatively 



the same as the situation described by Pontin et al. (2005) in the case of 
p — 1. The nature of the plasma flow in a plane of constant x with different 
values of p is shown in Fig. |2j Note that the qualitative structure - of a 
stagnation-point flow - is not affected a great deal by varying p. However 
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(a) (b) (c) 



Figure 2: Structure of the plasma flow across the spine and fan (black lines) in typical 
plane of constant x = 0, where the grayed area is the diffusion region, for (a) p = 2, (b) 
p = 0.9, (c) p = 0.5, for parameters rjo = jiQ = Bo = j = a = b = L = 1. 

the general trend is that as p tends to zero, the plasma flow across the fan 
plane becomes weaker. We will return to discuss this behaviour below. 

3.2 Reconnect ion rate 

It is generally accepted that magnetic reconnection plays a fundamental role 
in many types of explosive astrophysical phenomena, for example solar flares. 
Yet what determines the reconnection rate is still a major problem and this is 
an important aspect of any reconnection model. In general, the reconnection 
rate in 3D is defined by the maximal value of 



E\\ ds (20) 
along any field line threading a spatially localised diffusion region D (e.g. 



Schindler et al.[ 1988|). By symmetry, in this case 

En ds (21) 



C2 

where the curve C2 lies along the shown in figure [3} Since the fan 

is a flux surface, the integral may equally well be performed along the curve 
CI, the curve CI lying in the fan perpendicular to B, see Fig. [3} Now, since 
the curve CI lies outside D and therefore, along it v x B = — E, we can write 

F — — ( vxB-dl (22) 
Jci 
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Figure 3: The curves CI and C2 joining two points on the x-axis, where the grayed area 
is the diffusion region, the arrows indicate the direction of field lines, for a=l 



from which it is clear that this reconnection rate measures the rate at which 
flux is transported across the fan surface by the flow in the ideal region 



(Pontin et al. 2005 



From equation (|21|), we have 
F = 




< p < 1, 
p > 1. 



(23) 



where is given by Equation (19) as before. Note as a point of verification 
that this reduces to the expression found by Pontin et al. (2005 ) when p — 1. 



Here we consider the dependence of the reconnection rate on the parameter 
p in two distinct cases — see Fig. |4} First we set the parameter j to be a 
constant, j = 2j say, which results in a current which is dependent on p. 
We then go on to consider the case where we set j = jo(p + 1), so that the 
current (J = 2B j / (L(p + 1)) x) is independent of p. 

We will also consider the effect, in each of these cases, of taking different 
values for the parameter a, which controls the dimensions of the diffusion 
region. When a = 1, the diffusion region is symmetric for all p, having 
circular cross-section in any plane of constant z. However, as stated above, 
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5 



4 - 



3 - 

Rcco. rate 





Figure 4: Dependence of the reconnection rate on p, where the solid curve at a = 1.5, 
dash-dotted curve at a = 1, long dashed at a = 0.5, for (a) j = 2j , (b) j = jo(p + 1), and 
parameters tjo = /i = Bq = jo = L = 1. 
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the boundary of the diffusion region intersects each of the 3 coordinate axes 
at 

f x = ±a 1 /P, y = ±a, z = ±b^/ 2 , < p < 1, . , 

\x = ±a 1 / p , y = ±a, z = ±6(p+ 1 )/ 2 p, p > 1. ^ ' 

Thus the diffusion region becomes asymmetric in the x?/-plane when a ^ 1 
and p ^ 1. It will be seen later that this property is advantageous when 
comparing with the results of a numerical simulation. 



3.2.1 Reconnection rate as p — > oo 



In the limit p — » oo, we observe from (24) that the diffusion region becomes 
approximately symmetric (exactly symmetric if a = 1). The following all 
holds for all values of a. For the two choices of dependence for our parameter 



j stated above, evaluating Eq. (23) we find 



lim F 

p— >oo 



j=2jo 



lim F 

p— >oo 



4jo£ ?7o 



(25) 



(26) 



=j'o(p+i) Lfj, 

Consider first the case where the parameter j is chosen to be a constant, 
j = 2j (so that the current J = 4:B j /(L(p + l)//o) x). The magnetic field 
in this case is B — > (2B /L)(0, y, —z) as p — > oo. That is, the magnetic field 
approaches a 2D X-point structure with zero current, and so the result above 
(F — > 0) is as expected. 

By contrast, the reconnection rate approaches a constant finite value as 
p — > oo when j = jo(p + 1) (so that J = 2B j / Lfi x). In this case the 
magnetic field is B — > (2B /L)(0,y — jaz, —z). So the configuration is that 
of a 2D X-point with a uniform current (proportional to jo)- As the diffu- 
sion region has only a finite extent along the direction of the current (x), 
the reconnection rate is finite. Note that as expected it is proportional to 
the parameters i]o,Bq/L and jo, where t]q is the resistivity at the null, and 
2Boj /Lfi is the current modulus. 



3.2.2 Reconnection rate as p — > 

We now turn to the opposite limit; p — > 0. Note that our two parameter 
choices j = 2j and j = j (p + 1) clearly reduce to the same situation (with 
jo replaced by 2j ) as the limit is approached. Setting p = the magnetic 
field is B = (2Bq/L)(x, —joz, —z). We note that this field contains a neutral 
line in 3D (along y = 0) which is anti-parallel to the direction of current 
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flow — not a configuration associated with 2D reconnection. In fact the limit 



of Eq. (23) is not well defined for all choices of our parameters. Therefore 
we consider that p = is not a physically relevant parameter choice and 
consider only the limit p — > 0. 

As p — > 0, the magnetic field in the fan plane parallel to the current vector 
becomes strong, while the y-component becomes weak. Correspondingly, the 
flow across the fan surface becomes isolated to a small region near the fan, 
and weakens, see Fig. [2] Furthermore, in this case the diffusion region D is 
highly anti-symmetric. 



Evaluating Eq. (23) we find 




1 a > 1 

(taking j = 2 jo). For a < 1 the extent of D along the x-axis (direction 
of current flow) shrinks to zero. The result of the weak flow across the 
fan for small p is therefore that the reconnection rate also approaches zero 
when p — Y 0. By contrast, for a > 1 the boundaries of D stretch to infinity 



along x (see Eq. (24)). Correspondingly, for a > 1 the reconnection rate 
F — > oo. Although the flow is still very weak across the fan, the diffusion 
region now has much larger extent in the x-direction, and so although the 
flux reconnected per unit time per unit length in that direction decreases, 
the total flux reconnected increases. When a — 1, D is symmetric, because 
the boundary is at x = ±a, and the reconnection rate approaches a constant 
value as p — > 0. 

The results discussed above show that depending on our choice of pa- 
rameters there are various different ways in which the reconnection rate may 
depend on the asymmetry of the field (p). We now go on to perform simu- 
lations in the resistive MHD regime, in order to investigate which of these 
dependencies is relevant in a dynamically evolving plasma. 



4 Resistive MHD simulations 
4.1 Computational setup 

We now proceed to test the results of the mathematical model presented in 
the previous section by performing numerical simulations which solve the full 
set of resistive MHD equations. We solve the MHD equations in the following 
form 

E = -v x B + rjJ (28) 
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d 
dt 



J 

dB 

~dt 
dp 

dt 

(pv) 

de 
dt 



V x B 
-V x E 

-V ■ (pv) 

-V- (ppv + r) — VP + J x B 
-V ■ (ev) - PV ■ v + Q v i sc + Qj, 



(29) 
(30) 

(31) 

(32) 

(33) 



where v,B,E,rp J, p, r, P, e, Q„j SC , Qj are the velocity, magnetic field, electric 
field, resistivity, electric current, density, viscous stress, pressure, internal en- 
ergy, viscous dissipation and Joule dissipation, respectively. Here we provide 
a brief explanation of the method used for the numerical simulations. We run 



simulations that are similar to those described by Pontin et al. (2007). For 
more details on the numerical method, see Nordlund and Galsgaard (1997); 



Pontin and Galsgaard (2007). All simulations use numerical resolution of 
128 3 grid cells, a uniform resistivity model, and a so-called 'hyper viscos- 
ity' model. This is calculated using a combined 2nd- and 4th-order method, 
which effectively 'switches on' the viscosity only where grid-scale features 
develop in v, in order to maintain code stability. In this way the effect of vis- 



cosity is minimised, and we focus on the effect of the resistivity (see Nordlund 



and Galsgaard (1997)). 



We consider an isolated three dimensional null point within our compu- 
tational volume, which is driven from the boundary. We begin initially with 
a potential magnetic field 



B = -° — ^— (x,py, -(p + l)z) 
L p + 1 



(34) 



taking Po=L=l and ^=0.0007, constant, throughout. The computational 
domain has dimensions [x,y, z] = [— 3 . . . 3, — 3 . . . 3, — 0.5 . . . 0.5], with the 
magnetic field being line tied on all boundaries. In the kinematic model the 
spine and fan plane are not orthogonal, but in the our simulation at the out- 
set they are orthogonal, which means the plasma is in equilibrium. At t = 0, 
the spine of the null point lies in ^-direction, and the fan plane in the z = 
plane. A driving velocity is then assumed on the z-boundaries, which advects 
the spine footpoints in opposite directions on the opposite boundaries. The 
spine is driven until the resulting disturbance reaches the null, resulting in 
the formation of a current sheet as the magnetic field becomes stressed and 
distorted. The resulting configuration shares key properties with the config- 
uration considered in the kinematic model: the spine and fan have collapsed 
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(a) 



(b) 



(c) 



Figure 5: Structure of the magnetic field for the run with p = 2 (a) at t = 0, and (b) 
at the time of maximum current (t — 3.0), once the magnetic field has locally collapsed 
to form a current sheet. The black field lines are traced from around the spine for z > 0, 
while the grey field lines are traced from z < 0. (c) Grayscale showing |J| in the x = 
plane at the same time as (b). 

toward one another generating a current parallel to the fan surface, and fur- 
thermore a localised non-ideal region is present around the null. After some 
time the driving velocity is reduced back to zero. The explicit form taken for 
the driving velocity is defined by the streamfunction 

= ±0.5) = ±0.01 ((^H) 4 - l)\m (f ) cos 2 (f ) e -W\ 

(35) 

. Below we compare 
the kinematic solution and our simulation results. 



< t < 3.6 (for more details, see Pontin et al. (2007 



4.2 Current sheet 

In order to simplify the discussion we will initially explain the behaviour of 
the current at one value ofp(p = 2). We first examine the temporal evolution 
of current in the volume. In the beginning the spine and fan are orthogonal, 
but then the angle between them begins to change, reaching a minimum 
value once the current sheet forms. In other words the null collapses from a 
perpendicular X- type null point, with the angle between the X becoming 
greatly reduced, see Fig. [5j After the boundary driving ceases the current 
begins to decrease again, and the spine and fan relax back towards their 



initial perpendicular state, see Fig. 7(a) 
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We now discuss how the current sheet formation, as described above, 
depends on the value of p. Fig. [6] illustrates the dimensions of the current 
sheet for various values of p at the time when the current modulus is of 
maximum value. For the case investigated previously, p = 1, the sheet was 
found to be approximately of equal dimensions along x and y, the two coor- 
dinate directions associated with the fan surface. Looking at figure [6] we see 
a large difference between the geometry of the current sheet at p = 0.1 and 
at p = 10 at the maximum current. We find the current sheet at p = 0.1 is 
large, being very extended along the x-axis, that is, the direction along which 
J and the parallel electric field lie. However, this length decreases when p is 
increased. The results suggest that when the value of p approaches zero, the 
current sheet will grow indefinitely in the plane perpendicular to the shear, 
i.e. the direction of current flow through the null (x-direction). Note that 
with respect to the field strength in the fan plane, decreasing p corresponds 
to weaker magnetic field strength along the x-direction. Thus the extension 
of the current sheet could be attributed to the fact that the weak field re- 
gion extends in that direction and the magnetic field becomes less able to 
resist the collapse to form the current layer. That is, when the magnetic 
field parallel to the current becomes weaker there is less magnetic pressure 
associated with this 'guide field' component in the current sheet away from 
the null, and the current sheet is able to extend further away from the null. 



4.3 Maximum current attained 



Figure 7(a) illustrates the evolution of the current modulus maximum within 
the domain in time, for runs with different values of the parameter p. We 
notice in each case that the peak current grows sharply in time to reach a 



maximum value, and then decreases again, as discussed above. Figure 7(b) 
shows the maximum value in time of current plotted against p. We observe 
that there is a positive correlation between p and the maximum current - 
in other words, when p increases then the maximum value of the current 
that is attained also increases. Furthermore there is a negative correlation 
between the size of the current sheet and the value of p, see figure (|6]). 
Thus, although the current becomes more localised as p increases, it also 
becomes more intense. One important point to note is that of course the 
effective value of p will change during the simulations as the magnetic field 
is deformed. This can be confirmed be calculating the eigenvalues of VB at 
the null as the simulations proceed. We find that the relative change is small 
- of order 1% for p = 0.1 and order 10% for p = 10. Thus the ordering of the 
values of p that we selected for our simulations is preserved and the trends 
that are observed for the p-dependence will be unaffected. 
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(a) 



(b) 



(c) 



(d) 



(e) 



Figure 6: Isosurfaces of |J|, at 50% of the maximum at that time, showing dimensions 
of the current sheet for different p at the time of maximum |J|: (a) p = 0.1 (t = 3.8), (b) 
p = 0.5 (t = 3.5), (c) p = 1 (t = 3.3), (d) p = 2 (t = 3.0) and (e) p = 10 (i = 3.2). Plot 
dimensions are [x, y, z] = [±2.7, ±1, ±0.17]. 



4.4 Reconnection rate 

The nature of the plasma flow, and the resulting qualitative structure of 
the reconnection process, are found to be independent of the value of p. 
Specifically, we find plasma flow across both the spine line and fan plane of 
the null for all values of p. Figure [8] shows the plasma flow for two different 
values of p. Comparing with Figure |2j we see that the trend for the geometry 
of the flow is the same as in the kinematic solution. Specifically, for large p, 
the flow exhibits a relatively symmetric stagnation structure (in the x = 
plane). For smaller p the flow across the fan becomes confined to a narrower 
region, and comparatively weaker with respect to the flow across the spine. 

In this section we calculate the reconnection rate, i.e. the amount of 
flux transported across the fan surface, as before by integrating the electric 
field component parallel to the magnetic magnetic field (Eh). Similarly to 
above, by symmetry, the integral is performed along the field line lying along 
the x-axis, where since we are in the resistive MHD regime we have E» = 
77J • B/|B|. 
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Figure 7: (a) Evolution of the maximum value of |J| in time with different p; p = 0.01 
(dotted), p = 0.05 (long dashed), p = 0.1 (dash dot dot), p — 0.5 (dash dot), p = 1 
(dotted), p — 2 (dashed) and p = 10 (solid), (b) Peak spatial and temporal value of |J| 
for different p. 



In figure [9] we show the evolution of the reconnection rate in time for 
different values of p. Initially the rate clearly stay constant (zero) in time, 
i.e. during the early evolution, between t = and t — 1. Later, it starts to 
develop until it gains its maximum value, and then begins to decrease. This 
follows the same pattern as the evolution of the current, being indicative of 
the fact that the null point collapses to form the current sheet and recon- 
nection occurs, and then the system relaxes once the driving ceases. It is 
clear from Fig. [9] that the maximum reconnection rate attained increases as 
the value of p is decreased. It is worth emphasising here that although our 
intuition tells us that there is positive correlation between current and recon- 
nection rate, by contrast in this study we notice the inverse is true, i.e. when 
the peak current increases the reconnection rate decreases. This is because 
the diffusion region stretches when p tends to zero in the direction where 
the E\\ lies. Therefore the rate increases even though the current decreases, 



since the integrand in Eq. (20) is non-zero over a much larger portion of the 
x-axis. 

If we finally compare our results with those of the incompressible model 



of Craig and Fabling (1998), we find their results differ from ours in terms 



of the dependence of the peak current on p. In particular, they found that 
(in terms of our parameters) the maximum current decreases when p goes to 
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infinity. This may be down to the very different geometries of the current 
sheet in the two models (the current sheet in their incompressible model is 
planar and extends to infinity in all directions along the fan for all values 
of p). However, it is of interest to note that we actually find the same 
dependence of reconnection rate on p, i.e. as p decreases the reconnection 
rate increases (since in fact we find a negative correlation between J max and 
the reconnection rate as p is varied). 



5 Conclusions and Discussion 

In this paper we have investigated the effect of the symmetry of the magnetic 
field on magnetic reconnection at an isolated null point. We concentrate on 



the so-called spine-fan mode of 3D null point reconnection (Priest and Pontin 



2009). In a future paper we will go on to consider the 'torsional spine' and 
'torsional fan' modes, which involve a current flowing parallel to the spine of 
the null point. 

In the first part of the paper we discussed a steady solution of a subset of 
the resistive MHD equations, where the magnetic null point was defined by 
T5=B' (x,py — jz, —{p + l)z). This magnetic field has current aligned to the 
fan surface of the null point, and Pontin et al. ( |2005 ) investigated this situa- 



tion in the non-generic symmetric case p — 1 (repeated eigenvalues). In this 
work we use p as a parameter. By necessity, as the dynamics of the system 
are not included in this steady-state kinematic solution, a current is imposed, 
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Figure 9: (a) Reconnection rate at different values of p, where the dotted curve is the 
reconnection rate at p — 0.1, the long dashed curve at p = 0.5, the solid curve at p = 1, 
the dashed dot curve at p — 2 and the dashed dot dot curve at p = 10. (b) Variation of 
the maxiumum reconnection rate with the parameter p. 



which has the same orientation at the null as found in the simulations (the 
orientation of J at the null has been shown to be the crucial quantity in de- 



termining the topological structure of the reconnection process (Pontin et al. 



2004, 2005)). In order to have a localised diffusion region around the null, we 



artificially imposed a localised resistivity. We found the nature of the plasma 
flow, and the resulting qualitative structure of the reconnection process, to 
be the same as found in the symmetric case. Specifically, we found plasma 
flow across both the spine line and fan plane of the null for all values of p. 

We then described the results of a computational resistive MHD simula- 
tion in which we investigated the nature of the MHD evolution for different 
values of the parameter p (the ratio of the fan eigenvalues). Since in this 
case the full set of MHD quations was solved self-consistently, we began with 
an equilibrium potential magnetic null point (with J = 0). The system was 
then driven away from this equilibrium in such a way as to induce a local 
collapse of the null leading to current sheet formation and spine-fan magnetic 
reonnection. The resulting configuration shares key properties with the ana- 
lytical solution: the spine and fan are non-orthogonal with a current flowing 
parallel to the fan surface, and a localised diffusion region is focussed at the 
null. Also, in both cases the flow in the T/z-plane exhibits a stagnation-point 
structure. There is agreement between the model and the simulations, in 
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that for large p the stagnation structure is relatively symmetric, while for 
smaller p the flow across the fan becomes confined to a narrower region, and 
weaker compared with the flow across the spine. 

One of the major results that arises from the sequence of simulations is 
that both the peak intensity and the dimensions of current sheet are strongly 
dependent on the symmetry/asymmetry of the field in the fan surface, or in 
other words on the value of p. In terms of the sheet dimensions, the length 
along the direction of current flow at the null increases when p goes to zero, 
i.e. the diffusion region is stretched in the x-direction when p tends to zero. In 
the kinematic solution it was also possible by choosing the correct parameters 
to have the diffusion region dimensions have such a p-dependence. In order 
for our kinematic solution to be physically relevant, this implies that the 
parameter a in our solution should be chosen such that a > 1. Furthermore, 
as there is little difference in the size of the diffusion region in z for different 
p in the simulations, we should take b = 1 in our mathematical model. 

In addition to the current sheet at the null, we examined the reconnection 
rate in both cases. In order to compare the results, in light of the discussion 
above, we consider the parameter regime a > 1 in the kinematic solution. 
When a > 1 the reconnection rate — > oo as p — > 0. On the other hand, 
as p — > oo the reconnection rate approaches either zero or a constant finite 
value, depending on whether the current falls to zero or remains fixed, re- 
spectively, as p is increased (see Figure [4]). Turning to the simulations, as 
shown in Figure [9] the reconnection rate indeed becomes very large as p — > 0, 
in agreement with the kinematic model. In addition, as p becomes large 
the current at the null falls, and the reconnection rate appears to asymp- 
totically approach some small value, also in agreement with the kinematic 
model. Whether this value is finite or zero is not possible to tell within the 
restrictions of the present simulations. 

The results of both the mathematical model and simulations reveal that 
the symmetry /asymmetry of the magnetic field in the vicinity of a null can 
have a profound effect on the geometry of any associated reconnection region, 
and the rate at which the reconnection process proceeds. 
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